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摘 要 : 针对 骨架 颗粒 的 多 层 物 理 结构 问题 ， 运 用 孔道 网 络 方法 、 热 质 传递 原理 和 多 尺度 理论 等 知识 ， 建 立 了 颗粒 堆积 多 孔 介 
质 干燥 的 多 尺度 多 层 结构 传 热 传 质 模型 。 以 稻谷 堆 为 典型 代表 进行 热风 干燥 试验 验证 ， 模 拟 分 析 了 颗粒 内 部 不 同 组 织物 理 特 
性 等 因素 对 干燥 过 程 的 影响 。 结 果 表 明 : 建立 的 模型 可 有 效 模 拟 稻谷 堆 干 燥 过 程 ， 稳 谷 颗 粒 内 存在 较 大 的 水 分 梯度 。 胚 扩散 
系数 对 干燥 过 程 的 影响 十 分 显著 ， 其 次 是 壳 扩 散 系 数 ， 衣 扩散 系数 影响 最 小 ， 较 小 的 胚 扩散 系数 可 将 湿 分 有 效 地 “囚禁” 在 胚 
内 。 
关键 词 ， 颗粒 堆积 多 孔 介 质 ; 干燥 ; 多 层 结构 模型 ， 多 尺度 ;扩散 系数 
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Abstract: Aiming at the problem of multi-layer physical structure for the skeleton of porous media, a multi-scale 
and multi-layer structural model of heat and mass transfer processes for drying of grain packing porous media was 
established by applying the pore network method and multi-scale theory. An experimental study on rice drying was 
conducted in order to validate this model. The simulation and experimental results indicated that the established 
model could explain the mechanical properties of rice drying well. There was a higher moisture gradient inside the 
rice grain. The diffusion coefficient of rice embryo played an important role in the drying process, and whose effect 
on drying was larger than the diffusion coefficient of rice hull and chaff. The moisture was imprisoned effectively 
inside the rice when the diffusion coefficient of rice embryo was very small. 
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0 引 了 化 学 、 气 象 学 和 高 能 物理 等 科学 的 核心 理论 之 一 口 。 
孔道 网 络 干燥 理论 是 在 上 个 世纪 90 年 代 初 ，Daian 
颗粒 堆积 多 孔 介质 干燥 涉及 农业 、 食 品 、 轻 “和 Saliba 以 及 Nowicki 等 人 率先 将 渗流 理论 中 常用 的 
工 、 化 工 、 建 筑 、 煤 炭 等 众多 领域 。 目 前 ， 对 农 产 ”孔道 网 络 方法 引入 到 多 孔 介质 干燥 理论 中 ， 后 经 Prat 
品 颗 粒 堆积 多 孔 介 质 干 燥 的 机 理 及 模拟 研究 ， 传 统 。 和 Yortsos 等 人 的 进一步 研究 完善 而 逐渐 发 展 起 来 的 
的 方法 主要 是 连续 介质 法 ， 即 将 这 类 物料 视 为 虚拟 。 ”1。 这 一 理论 认为 多 孔 介 质 的 湿 分 迁移 在 孔隙 中 具 
的 连续 湿 固 体 ， 不 考虑 物料 堆 中 颗粒 骨架 与 孔隙 的 。” 有 很 大 优势 ， 用 一 系列 规则 排列 的 节点 及 与 其 相互 连 
差别 ， 也 不 考虑 颗粒 内 部 组 织 结构 的 不 同 ， 直 接 运 接 的 孔道 来 重 现 多 和 孔 介质 和 孔 际 空 间 的 几何 及 拓扑 结 
用 现 有 热 质 传递 等 理论 进行 建 模 求解 中 。 而 事实 。 构 ， 直 接 在 孔道 等 级 上 研究 干燥 过 程 中 多 孔 介质 的 热 
上 上 ， 农 产品 颗粒 堆积 多 孔 介 质 结构 复杂 ， 在 大 部 分 。 质 传递 问题 。 
情况 下 更 接近 于 离散 介质 ， 且 其 骨架 颗粒 内 部 是 一 本 文 拟 将 多 尺度 理论 的 基本 思想 引入 颗粒 堆积 
个 多 层 物理 结构 ， 包 括 种 皮 《〈 壳 ) 、 种 衣 和 胚 体 等 ”” 孔 介 质 干 燥 研 究 领域 ， 综 合 运用 孔道 网 络 方法 、 热 质 
组 织 。 显 然 ， 忽 略 这 些 因素 对 于 燥 过 程 的 影响 是 不 传递 原理 等 交叉 学 科 的 知识 ， 建 立 其 干燥 过 程 的 多 尺 
符合 实际 情况 的 。 度 多 层 结构 传 热 传 质 模型 。 并 以 稻谷 堆 为 颗粒 堆积 多 
多 尺度 理论 是 研究 各 种 不 同 空间 尺度 或 时 间 尺 ”和 孔 介 质 的 典型 代表 进行 热风 干燥 试验 验证 ， 模 拟 分 析 
度 相互 耦合 现象 的 科学 。 经 过 近 三 十 年 的 发 展 ， 多 ”颗粒 内 部 不 同 组 织物 理 特性 等 因素 对 干燥 过 程 的 影 
尺度 理论 已 是 流体 动力 学 、 材 料 科 学 、 环 境 科 学 、 响 。 
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1 二 维 多 尺 度 多 层 结构 传 热 传 质 模型 


1.1 物理 模型 

农产品 颗粒 堆积 多 孔 介 质 通风 干燥 过 程 主要 
是 指 热风 在 风 压 力 的 推动 下 ， 从 物料 的 某 个 或 多 个 
方向 进入 其 孔 隐 体系 ， 与 骨架 颗粒 进行 热 质 交 换 后 
从 另 一 个 或 多 个 方向 排出 的 过 程 ， 如 图 la 所 示 。 
显然 ， 在 热风 干燥 过 程 中 孔隙 是 热 质 传递 的 优势 通 
道 。 为 区 分 颗粒 骨架 与 孔隙 的 差别 ， 本 文采 用 文献 
(I 中 的 方法 构建 孔道 网 络 物理 模型 如 图 lb 所 示 。 
模型 主要 由 “ 孔 - 喉 - 孔 ”网 络 和 骨架 两 部 分 组 成 ， 孔 
为 球体 ， 喉 为 圆柱 体 ， 孔 和 喉 径 服从 一 定 的 概率 分 
布 。 一 个 二 维 孔 道 网 络 物理 模型 可 由 分 布 律 、 两 节 
点 间距 离 、 模 型 规模 数 、 配 位 数 等 参数 来 描述 
9 。 其 次 ， 为 考查 骨架 颗粒 内 部 壳 、 衣 、 胚 等 组 
织 结构 不 同 对 干燥 过 程 的 影响 ， 本 文 将 研究 对 象 划 
分 为 颗粒 尺度 与 干燥 器 尺度 两 个 不 同 级 别 来 处 理 。 
在 颗粒 尺度 上 ， 每 一 个 骨架 颗粒 都 具有 种 皮 
《 壳 ) 、 种 衣 和 胚 体 三 层 不 同 的 物理 组 织 结构 ， 如 
图 lc 所 示 。 
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(a) 热风 干燥 (b) 孔道 网 络 “(Cc) 颗粒 尺度 物理 结构 
1 多 尺度 多 层 结构 模型 
Fig.1 Model of multi-scale and multi-layer 

1.2 数学 模型 

为 简化 数学 模型 ， 针 对 上 述 颗 粒 堆 积 多 孔 介 质 
热风 干燥 过 程 ， 我 们 作 如 下 基本 假设 : (1) 在 颗 
粒 尺 度 上 假设 骨架 颗粒 为 具有 多 层 物 理 结构 的 椭 球 
体 ， 不 收缩 、 无 变形 ; 假设 颗粒 内 部 湿 分 传递 主要 
依靠 液态 或 汽 态 扩散 ， 假 设 骨架 颗粒 内 的 水 分 相 变 
所 需 耗 热 由 骨架 颗粒 本 和 喘 提供 。 (2) 在 干燥 器 尺 
度 上 以 对 流 换 热 为 主 ， 忽 略 热传导 、 辐 射 传 热 ， 对 
流传 质 为 主 ， 忽 略 扩散 传 质 ， 不 考虑 有 液 相 参与 的 
各 种 传 质 ;假设 气相 为 不 可 压缩 的 理想 气体 。 
1.2.1 颗粒 尺度 

对 于 孔道 网 络 物理 模型 中 在 颗粒 尺度 上 的 每 个 
骨架 颗粒 (图 lc) ， 可 采用 雷 可 夫 理 论 " 描述 其 
内 部 的 热 质 传递 过 程 ， 由 上 述 基本 假设 ， 经 简化 推 
导 可 得 在 圆柱 坐标 系 下 的 热 质 传递 方程 组 为 
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式 中 : 7 为 颗粒 内 部 某 时 刻 某 处 的 温度 ，K ;4 为 颗 
粒 热 导 率 ，W/m:KK; ps 为 颗粒 密度 ，kg/m ; C。 为 颗 
粒 比 热 容 ，J/kg:K; M 为 颗粒 内 部 某 时 刻 某 处 的 湿 
含量 (db.) ;也 为 颗粒 内 水 分 扩散 系数 ，mYs; t 为 
时 间 变 量 ，s; r、2Z 为 圆柱 坐标 变量 ，m。 
1.2.2 干燥 器 尺度 
如 图 1b， 设 在 某 一 时 刻 nA t， 了 和 孔道 网 络 物理 模 
型 中 某 孔 节点 (i, 妃 内 的 燕 汽 密度 为 p,"y ， 则 应 用 质 
量 、 能 量 守 恒定 律 与 达 西 定律 可 得 在 干燥 器 尺度 上 的 
传 质 、 传 热 及 动量 传递 控制 方程 组 为 
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式 中 : pi 为 某 一 时 刻 某 孔 节 点 仿 刀 内 的 蒸汽 密 
度 ，kg/m’; Rg 为 节点 也 半径，m; v= g; /人 (Mm)， 
为 喉 道 中 风速 , m/s，; i 为 喉 半 径 ，m; 4s 为 颗粒 的 
表面 积 ，m2; 到 为 对 流传 质 系 数 ，m/s;， pe 为 颗粒 表 
面 的 蒸汽 密度 ，kg/m ;At 为 时 间 步 长 ，s，7 "为 革 
一 时 刻 某 孔 节 点 (i, 站 内 的 气相 温度 ，K; p 为 气相 密 
度 ，kg/m ; C 为 气相 比热容 ，J/Kg.KK; hh 为 对 流 换 
热 系数 ，W/(m-K); gi 为 喉 中 的 气相 体积 流 率 ， 
m/s; 4; 为 喉 截 面积 ，m A* 为 达 西 渗透 率 ， m’; 从 
为 气相 粘度 ，Pa.s; ] 为 喉 长度 ，m; PP 为 孔 节 点 风 
压 ，Pa; pw 为 纯 蒸汽 密度 ， kg/m ; 
Sk = horke(Pr ot 一 Pr 1)14， 为 颗粒 向 某 孔 释放 的 
蒸汽 速率 ，kg /s。 
1.2.3 跨 尺 度 热 质 交 换 控 制 方程 

上 述 两 个 尺度 方程 之 间 的 热 质 传递 信息 彼此 分 
离 ， 没 有 联系 ， 故 需 将 两 者 的 信息 联系 起 来 ， 补 充 一 
定 的 跨 尺 度 热 质 交 换 平 衡 方程 方 可 求解 。 

在 孔道 网 络 物理 模型 的 每 个 骨架 颗粒 表面 上 ， 存 
在 着 跨 尺 度 的 热 质 交换 ， 即 该 颗粒 内 水 分 的 减少 量 等 
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于 气 化 流入 颗粒 周围 孔隙 的 蒸汽 量 ， 同 时 ， 该 颗粒 
内 热 烩 变化 量 等 于 颗粒 从 周围 气相 中 对 流 换 热 得 到 
的 热量 减 去 气 化 颗粒 内 水 分 所 耗 的 热量 。 于 是 应 用 
人 


质量 、 能 量 守 恒定 律 可 得 跨 尺 度 的 热 质 交 换 平 衡 方 
程 组 为 
4xRp,AX 之 4, 
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式 中 : 0 为 颗粒 绝 干 密度 ，kg/mm we 为 水 分 气 化 
潜 热 系数 ，J/kg; 人 站 为 某 颗 粒 在 t 时 间 内 的 含水 
率 差 ， Kg/ kg。 
2 验证 性 干燥 试验 

以 稻谷 扒 为 颗粒 堆积 多 孔 介质 的 典型 代表 进行 
热风 干燥 试验 验证 。 首 先 ， 将 稻谷 除 杂 洗 净 并 调 质 


获得 均匀 一 致 的 湿 含量 ， 然 后 ， 将 稻谷 置 于 干燥 斌 
验 仓 中 进行 热风 干燥 ， 同 时 测 得 稻谷 堆 主要 结构 特 


征 参 数 。 干 燥 仓 为 宽 200mm， 高 200mm， 厚 
100mm 的 长 方 腔 体 ， 仓 底部 中 心 处 开 一 直角 等 腰 


三 角形 进 风 口 ， 顶 部 敞开 ， 与 外 界 环境 空气 接触 ; 
在 仓 壁 正面 开 有 测量 孔 。 


日 干 || 
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中 ， 干 燥 模 型 的 初始 、 边 
致 。 


结果 分 析 与 讨论 
4.1 干燥 仓 内 稻谷 平均 湿 含量 变化 曲线 

从 图 2 可 以 看 出 ， 仓 内 稻谷 颗粒 的 平均 湿 含 量 模 
拟 与 实验 变化 曲线 基本 吻合 ， 其 最 大 相对 误差 约 为 
7.6%。 这 从 宏观 上 表明 ， 本 文 建立 的 模型 有 效 。 
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图 2 仓 内 稻谷 颗粒 的 平均 湿 含 量 曲线 模拟 与 实验 结果 的 比较 
Fig.2 Drying curves comparison of simulation and experiment 
4.2 颗粒 内 部 不 同 组 织物 理 特性 对 干燥 过 程 影 响 
4.2.1 壳 扩散 系数 

将 壳 扩 散 系 数 设 定 为 5 种 不 同 的 值 ， 其 余 参数 不 
变 ， 进 行 颗粒 堆积 多 孔 介质 干燥 模 拟 ， 得 到 不 同 壳 扩 
散 系 数 下 的 颗粒 干燥 曲线 如 图 3 所 示 。 从 图 中 可 以 看 
出 ， 颗 粒 壳 扩散 系数 对 干燥 过 程 有 一 定 影响 。 随 着 碗 
扩散 系数 的 减 小 ， 干 燥 速 度 逐 渐变 慢 ， 达 到 相同 湿 含 


干燥 试验 开始 时 ， 仓 进 风 口 热风 温度 
42.0C， 热 风 相 对 湿度 6.5%， 热 风 静 压 56.0Pa; 
环境 温度 20.1*C; 稻谷 颗粒 初始 温度 18.0*C， 稳 
谷 颗粒 初始 湿 含量 31.2% (db.)。 于 燥 试 验 过 程 
， 每 隔 一 段 时 间 拔 开 测量 孔 橡胶 塞 ， 用 红外 测 温 
仪 测量 各 点 稻谷 颗粒 的 温度 ， 并 从 测量 孔 中 取出 适 
量 稻谷 样品 ， 采 用 烘 干 法 测 得 各 样品 的 湿 含量 。 仓 
内 稻谷 堆 中 孔隙 气相 的 温 、 湿 度 由 埋 入 各 测量 点 的 
温 、 湿 度 传感器 测量 


3 计算 机 模拟 


基于 上 述 干 燥 模 型 ， 本 文采 用 VC++ 与 
MATLAB 语言 联合 编程 开发 模拟 程序 。 干 燥 器 尺 
度 控制 方程 和 跨 尺 度 热 质 交换 平衡 方程 已 是 离散 型 
的 差分 方程 ， 直 接 采 用 有 限 差分 的 分 离 式 解法 求 
解 ;而 颗粒 尺度 方程 为 偏 微分 形式 ， 需 先进 行 网 格 
划分 然后 对 其 有 限 差 分 离散 方 可 求解 。 稻 谷 颗粒 为 
se 为 提高 计算 速度 ， 只 
需 取 其 对 称 面 的 四 分 行 分 析 。 模 拟 计算 过 程 


本 


量 所 需 干 燥 时 间 越 长 。 这 主要 是 因为 帝 扩 散 系数 越 
小 ， 颗 粒 内 湿 分 逸 出 壳 表 面 越 困 难 ， 王 燥 也 就 越 慢 。 
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3 不 同 壳 扩散 系数 的 颗粒 干燥 曲线 


Fig.3 Drying curves of simulation for different hull diffusion 


coefficient 

4.2.2 衣 扩 散 系 数 

将 衣 扩 散 系 数 设 定 为 5 种 不 同 的 值 ， 其 余 参数 不 
变 ， 进 行 颗粒 堆积 多 孔 介 质 干燥 模拟 ， 得 到 不 同 衣 扩 
散 系 数 下 的 颗粒 干燥 曲线 如 图 4 所 示 。 从 图 中 可 以 看 
出 ， 颗 粒 衣 扩散 系数 对 干燥 过 程 的 影响 很 小 。 随 着 衣 
扩散 系数 的 减 小 ， 干 燥 速 度 虽 有 所 变 慢 ， 但 变 慢 的 幅 
度 很 小 ， 可 以 忽略 不 计 。 这 主要 是 因为 颗粒 种 衣 的 组 
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织物 理 厚 度 相 对 这、 及 来 说 注 得 多 ， 其 扩散 系数 的 
变化 对 颗粒 内 部 湿 分 迁移 所 起 的 作用 有 限 。 
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4.2.3 胚 扩散 系数 
将 胚 扩散 系数 设 定 为 5 种 不 同 的 值 ， 其 余 参 数 
不 变 ， 进 行 颗粒 堆积 多 孔 介 质 干燥 模拟 ， 得 到 不 后 
及 扩 散 系 数 下 的 颗粒 干燥 曲线 如 图 5 所 示 ， 湿 含量 
分 布 (t=300 min 时 ) 如 图 6 所 示 。 
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diffusion coefficient 

由 图 5 可 以 看 出 ， 颗 粒 胚 扩散 系数 对 干燥 过 程 
的 影响 巨大 。 随 着 胚 扩散 系数 的 减 小 ， 干 燥 速 度 明 
显 变 慢 ， 特 别 是 干燥 的 中 后 期 ， 达 到 相同 湿 含量 所 
需 干 燥 时 间 显 著 增 加 。 这 一 方面 是 因为 种 豚 在 颗粒 
组 织 中 占有 较 大 比例 ， 另 一 方面 ， 从 图 6 也 可 以 看 
出 ， 胚 扩散 系数 越 小 ， 颗 粒 内 湿 分 从 胚 内 迁移 到 壳 
表面 越 困 难 ， 湿 分 仿佛 被 “办 禁 ?在 胚 内 ， 如 当 胚 扩 
散 系 数 为 2.98x101 m2s 时 ， 其 颗粒 中 心 区域 的 湿 
含量 几乎 没有 减 小 。 因 此 ， 干 燥 也 就 越 慢 。 
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5 结论 


本 文 建立 的 多 尺度 多 层 结构 传 热 传 质 模型 可 有 效 
模拟 颗粒 堆积 多 孔 介质 干燥 过 程 。 模 拟 得 到 的 仓 内 稳 
谷 干燥 曲线 反映 了 其 干燥 过 程 的 实际 情况 。 颗 粒 内 部 
不 同 组 织物 理 特 性 对 颗粒 堆积 多 孔 介质 干燥 有 较 大 的 
影响 ， 胚 扩散 系数 的 影响 十 分 显著 ， 其 次 是 壳 扩 散 系 
数 ， 衣 扩散 系数 影响 最 小 。 在 干燥 过 程 中 ， 随 着 胚 扩 
散 系 数 的 减 小 ， 干 燥 速 度 明 显 变 慢 ， 特 别 是 干燥 的 中 
后 期 ， 达 到 相同 湿 含 量 所 需 干燥 时 间 显著 增加 ; 较 小 
的 胚 扩散 系数 可 将 湿 分 “办 禁 ”在 胚 内 ， 使 湿 分 较 难 从 
颗粒 内 迁移 到 壳 表 面 。 
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